Probing topological order with Renyi entropy 
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We present an analytical study of the quantum phase transition between the topologically ordered toric-code- 
model ground state and the disordered spin-polarized state. The phase transition is induced by applying an 
external magnetic field, and the variation in topological order is detected via two non-local quantities: the 
Wilson loop and the topological Renyi entropy of order 2. By exploiting an equivalence with the transverse- 
field Ising model and considering two different variants of the problem, we investigate the field dependence of 
these quantities by means of an exact treatment in the exactly solvable variant and complementary perturbation 
theories around the limits of zero and infinite fields in both variants. We find strong evidence that the phase 
transition point between topological order and disorder is marked by a discontinuity in the topological Renyi 
entropy and that the two phases around the phase transition point are characterized by its different constant 
values. Our results therefore indicate that the topological Renyi entropy is a proper topological invariant: its 
allowed values are discrete and can be used to distinguish between different phases of matter. 
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I. INTRODUCTION 

There are states in quantum many-body physics that can- 
not be described in terms of local order parameters and the 
Landau paradigm of spontaneous symmetry breaking. These 
states exhibit a subtler kind of order called topological order 
01]. Topologically ordered states include fractional quantum 
Hall liquids l2l and quantum spin liquids [3], which are at the 
forefront of research in condensed matter theory. Moreover, 
such states are of great interest in the field of quantum com- 
putation because one can encode quantum information in the 
topological degrees of freedom and this way of encoding is 
intrinsically robust against decoherence HIS]- 

Since topologically ordered states cannot be characterized 
by local order parameters, there has been an intense effort 
to find non-local quantities that can detect topological order 
in a wavefunction. A series of papers suggested that topo- 
logical order can be detected through a component of quan- 
tum entanglement that contains a topological constraint. This 
constraint manifests itself as a universal negative correction 
to the boundary law for the entanglement entropy: the so- 
called topological entropy Recent works have showed 
that this component of quantum entanglement is indeed long 
ranged and so it cannot be destroyed by time evolution with a 
local Hamiltonian. Equivalently, the components correspond- 
ing to the long-range entanglement and the usual short-range 
entanglement are adiabatically disconnected fioll . 

These results suggest that the topological entropy is to some 
extent a non-local order parameter for topologically ordered 
phases. To make a more precise statement about the extent of 
its applicability, one needs to investigate its robustness against 
perturbations. If the topological entropy is to detect topolog- 
ically ordered phases, it needs to be non-zero within all such 
phases. In other words, it should only vanish at quantum phase 
transitions to disordered phases. If it is to distinguish differ- 
ent topologically ordered phases from each other, it needs to 
be constant within each phase. In other words, it should only 
change at quantum phase transitions Hill . 

Recent numerical studies on small systems have found ev- 



idence that the topological entropy takes discrete values ^ 
and only changes at quantum phase transitions lfl2ll . On the 
other hand, analytic corrections to the topological entropy are 
extremely hard to obtain because one needs to consider the en- 
tanglement entropy of a many -body wavefunction Il3ll . Such 
a treatment for the topological entropy in the case of a finite 
correlation length can be found in Ref. Hill . Remarkably, 
it has been recently shown for an exactly solvable two-phase 
system that the topological entropy is constant within the en- 
tire topologically ordered phase 111 511 . 

Given the difficulties, it is important to find other entropic 
quantities that possess a topological component capable of de- 
tecting topological order. One potential candidate is the Renyi 
entropy of order a, which is a generalization of the usual 
(von Neumann) entanglement entropy. It is important that the 
Renyi entropy coincides with the entanglement entropy in the 
special case of a = 1. It has also been shown that the ground- 
state Renyi entropies of different order a all contain the same 
topological component at the fixed points of non-chiral phases 
J8fl. Such phases typically appear in string-net models lfl6ll 
and quantum double models 111711 . 

In this paper, we consider the Renyi entropy of order 2 
and argue that it is a good probe of topological order because 
its topological component can only change at quantum phase 
transitions. In particular, we apply the concept of the topolog- 
ical Renyi entropy to the toric-code model (TCM) |4(j in the 
presence of an external magnetic field. This model is, to para- 
phrase Goldenfeld lll8tl . the Drosophila of topological order. 
Although it is a simple toy model, it contains all the elements 
that make topological order interesting: there is no local order 
parameter, there is a topology-dependent ground-state degen- 
eracy that is robust against local perturbations, and there are 
excitations with anyonic particle statistics. Indeed, the TCM 
is another beautiful example of the crucial role played by toy 
models in statistical mechanics. 

To show that the topological Renyi entropy is a good probe 
of topological order, we demonstrate that the disordered and 
the topologically ordered phases of the TCM with external 
magnetic field are characterized by its distinct values. We also 
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study the Wilson loop as a probe of both topological order and 
gauge structure. Concentrating on two different variants of the 
problem, we establish an exact treatment in the computation- 
ally simpler (exactly solvable) variant and supplement it with 
perturbation theories in both variants. The results obtained 
with the two methods for the two probing quantities in the 
two variants are highly consistent with each other. 

H. GENERAL FORMALISM 

We consider the TCM with an external magnetic field in 
the +z direction. The system is an N x N square lattice with 
periodic boundary conditions, and 2N 2 spins are located at 
the edges of the lattice 0]. In general, the spins on the hori- 
zontal (h) and the vertical (v) edges experience different mag- 
netic fields: A on the horizontal and k\ on the vertical edges 
(k > 0). The Hamiltonian of the system takes the form 

H = -^i s -^4-A^af- K A^^, 

s p ieh i£v 

A s = l[af, B p = \{at, (1) 

where the indices s and p refer to stars and plaquettes on the 
lattice containing four spins each. For an illustration of this, 
see Fig. Q] Note that the four sums in Eq. ([TJ all contain 
TV 2 terms, and that only TV 2 — 1 star (plaquette) operators are 
independent because Y[ s A s = JT B p = 1. 




FIG. 1: (Color online) Illustration of the square lattice with the phys- 
ical spins located at the horizontal (black circles) and the vertical 
(white circles) edges. Examples of a star (red cross labeled s) and a 
plaquette (blue square labeled p) are included. 

The TCM with zero external field (A = 0) is exactly solv- 
able because the stars A s and the plaquettes B p all commute 
with each other. The ground state is fourfold degenerate: there 
are four linearly independent states with A s = (A s ) = +1 
and Bp = (B p ) = +1 for all s and p. These degenerate 
ground states are distinguished by the topological quantum 
numbers Z\ = ±1 and Z2 = ±1, which are expectation 
values for products of <xf operators along horizontal and ver- 
tical strings going round the lattice. The ground state with 
Z\ = Z-2 = +1 can be written as 

\0)=M]j(i + A s )\it), (2) 



where Af = l/y / 2 N2 + 1 is a normalization constant, and | ff ) 
denotes the completely polarized state with all spins pointing 
in the +z direction (af = +1 for all i). 

The TCM with finite external field (A > 0) is not exactly 
solvable because the magnetic fields af do not commute with 
the stars A s . On the other hand, they commute with the pla- 
quettes Bp and the topological operators Z12, therefore there 
are four independent lowest-energy sectors with Z\^ = ±1 
and Bp = +1 (Vp). In the rest of the paper, we consider the 
lowest-energy eigenstate |f2(A)) within the Z\ = Z2 = +1 
sector. This state becomes | ) in the limit of A — > 00 and 
|0) in the limit of A = 0. Between the two limits, numerical 
studies reveal a quant um p hase transition at a critical mag- 
netic field A = Ac lfl2lll9ll . Since |f2(A)) i s a ground state at 
the fixed points of both limiting phases, the adiabatic theorem 
guarantees that it is the unique ground state in the disordered 
phase at A > Ac and one of the four degenerate ground states 
in the topologically ordered phase at A < Ac- 

If we only consider the states with Z\ = Z- 2 = +1 and 
Bp = +1 (Vp), the dimension of the effective Hilbert space 
is reduced from 2 2N to 2 N _1 . The states within this re- 
duced Hilbert space can be written as superpositions of loop 
configurations on the dual lattice: in each loop configuration, 
the spins on the loops have af = — 1 and the remaining spins 
have af = +1. This implies that the reduced model is equiv- 
alent to a Z2 lattice gauge theory, and the phase transition at 
A = Ac corresponds to a confinement-deconfinement tran- 
sition IU2l [l9tl . Furthermore, since each loop configuration 
can be characterized by the values of the iV 2 — 1 independent 
stars A s = ±1, it is convenient to introduce a corresponding 
representation in which quasi-spins A s are located at the stars 
Il20ll . This quasi-spin representation is particularly useful in 
the A <C 1 limit because |f2(A)) is then close to |0) which is a 
product state with A s = +1 (Vs). Up to an irrelevant additive 
constant, the Hamiltonian in Eq. (HJ becomes 

s {s,s')eh (s.s')Gv 

where (s, s') means that the summation is over horizontal and 
vertical edges between nearest-neighbor stars s and s'. Note 
that A z s = A s measures and A x s switches the quantum number 
A s , therefore the quasi-spin operators A z s and A x s satisfy the 
standard spin commutation relations. In the quasi-spin rep- 
resentation of Eq. (O, the TCM with external magnetic field 
is equivalent to a two-dimensional (2D) transverse-field Ising 
model (TFIM) in which the coupling strengths on the horizon- 
tal and the vertical edges are different in general. 

m. MEASURES OF TOPOLOGICAL ORDER 

We aim to describe how the topological order in the ground 
state |0(A)) changes as a function of A between the topo- 
logically ordered limit at A = and the disordered limit at 
A — > 00. To quantify topological order in an analytically 
tractable manner, we consider two measures: the Wilson loop 
and the topological Renyi entropy. 
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A. General properties 

The Wilson loop for a region R on the dual lattice is defined 
as the expectation value Wr of the operator 

w R = n *f = n m> w 

iG9fl sG-R 

where <9i? denotes the boundary of i?. If the region R is 
macroscopic with linear dimension D ^> 1, the Wilson loop 
follows a perimeter law Wr oc cxp(— /?£)) in the presence of 
topological order and an area law Wr oc exp(— (3D 2 ) in the 
absence of topological order lfl2ll . In this paper, we assume 
that the region R is a D x D square (see Fig. |2j. 
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FIG. 2: (Color online) Illustration of the Wilson loop for a square 
region with D = 3. The region R contains D 2 stars (red crosses) and 
the boundary OR (dashed line) contains 4D spins (yellow circles). 
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FIG. 3: Illustration of the subsystems A and B in the four cases 
(m) that are used to calculate the topological Renyi entropy. Each 
subsystem A has extension D and thickness d with J) > d> 1. 

The topological Renyi entropy is based on the paradigm 
of quantum entanglement. The Renyi entropy of order a be- 
tween two complementary subsystems A and B = A reads 

St B = log 2 Tr \pl\ = — !— log 2 Tr [p%] , (5) 

1 — a 1 — a 

where and ps are the reduced density operators for A and 
B. The topological contribution to the Renyi entropy can be 
extracted by taking a suitable linear combination of Renyi en- 
tropies that are calculated for different choices of the subsys- 
tems A and B 0] • In fact, the standard definition for the topo- 
logical Renyi entropy of order a is 

S£ = -SW + SV> + SW-S<?>, (6) 



where Sa = S^ B in the four cases (to) of partitioning the 
system shown in Fig. [3] The characteristic linear dimensions 
are the extension D and the thickness d of the subsystem A in 
all cases. To obtain a meaningful topological measure, these 
dimensions need to be macroscopic (D > d ^> 1). 

The topological Renyi entropy is non-zero if and only 
if the given state exhibits topological order [0 Ht]. For the 
TCM with external magnetic field, £f£ = for the disordered 
ground state | ft ) and = 2 for the topologically ordered 
ground state |0). In this paper, we demonstrate that S J detects 
the presence of topological order in the entire topologically or- 
dered phase at A < Ac ■ Note though that is independent of 
a at the fixed point of a generic non-chiral phase JH . Since the 
fixed points of different topologically ordered phases do not 
necessarily have distinct values of 5j, the topological Renyi 
entropy is unable to provide a complete characterization of a 
topologically ordered phase. 



B. Z2 lattice gauge theory 

Since the TCM is perturbed with external fields erf that 
commute with the plaquettes B p and the topological operators 
Zi t 2, the ground state |f2(A)) belongs to the lowest-energy 
sector with Z\ = Z2 = +1 and B p = +1 (Vp) for all values 
of A. If we only consider this sector, the gauge structure con- 
straint B p \ty) = 1^) is enforced on all states therefore 
the model is equivalent to a Z2 lattice gauge theory. 

An arbitrary state within the gauge theory can be ex- 
pressed as a superposition of loop configurations. Each con- 
figuration is a finite set of closed loops on the dual lattice, and 
the spins on the loops are flipped with respect to the remain- 
ing ones. These properties motivate us to introduce a modified 
definition for the topological Renyi entropy in which the sub- 
system A is substituted by the boundary DA between A and B 
in each case (m) of partitioning the system. This means that 
= Sf a ^a = gdA in Eq Q Formally, we define C as 
the set of star operators acting on both subsystems A and B, 
and dA as the set of spins that are only acted upon by stars in 
C. For an illustration of this, see Fig. [4] The boundary dA is 
always a finite set of closed loops on the real lattice: the num- 
ber of loops is n = 2 in the cases ( 1 ) and (4) , while it is n = 1 
in the cases (2) and (3). Since a loop on the real lattice and a 
loop on the dual lattice can only intersect at an even number of 
points, there are an even number of spins flipped on each loop 
of dA. This topological constraint ensures that the modified 
has similar properties to the standard one. For example, it 
is still true that 5 J = for | ft ) and S% = 2 for |0). 

The calculations in the rest of the paper are immensely sim- 
plified by using the modified definition for S^- Since the 
group generated by the star operators acting exclusively on 
the boundary subsystem dA only contains the identity, the re- 
duced density matrix pqa is diagonal in the basis of the physi- 
cal spins erf. Each diagonal element (paA)ss gives the proba- 
bility that 1^) realizes a given spin configuration {S| = ±1} 
in dA. Equivalently, if we choose a random loop configura- 
tion according to the probability distribution given by the state 
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FIG. 4: (Color online) Illustration of the subsystems in case (1) with 
dimensions D — 6 and d — 2. Spins are either in subsystem A 
(black circles) or in subsystem B (white circles). Stars in the set C 
are marked by red crosses, and spins in the subsystem dA are marked 
by blue rectangles. The boundary contains n = 2 closed loops on 
the real lattice with a combined length L = 32. 



\i$f), the probability of the spin configuration {E| } in dA is 
P[{Sf}] = (p3a)se- If we then choose two random loop 
configurations according to the same distribution, the proba- 
bility of them having the same spin configuration in dA is 



Tr 



[pIa] 



(7) 



This result motivates us to consider the topological Renyi en- 
tropy of order 2. In terms of the probabilities V^ m ' in the four 
cases (to) of partitioning the system, this quantity takes the 

form S% = k>g a [7>W 4 VP (2) P (3) ]. 

We can now develop an intuitive understanding of the phase 
transition by considering the two limiting cases. In the topo- 
logically ordered ground state at A < 1, the spin loops are 
deconfined and all possible loop configurations are equally 
probable. This means that the allowed spin configurations in 
the subsystem dA also share the same probability: the inverse 
number of allowed spin configurations. It is important that the 
number of boundary loops is n = 2 in the cases (1) and (4), 
while it is n = 1 in the cases (2) and (3). The cases (1) and 
(4) are therefore more constrained and have less allowed spin 
configurations in dA. This implies V {1) , V (i) > p( 2 ),p( 3 ) 
and S-f > 0. More precisely, since the constraint on each 
boundary loop reduces the number of allowed spin configura- 
tions by a factor of 2, the topological Renyi entropy is given 



,.(!) 



,(2) 



,.(3) 



.(*) 



2. In the disordered 
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ground state at A ^> 1, the spin loops are confined and only 
the loop configurations with small spin loops have significant 
probabilities. On the other hand, the small spin loops in these 
loop configurations correspond to local disturbances (nearby 
spin flips) in the spin configurations of the boundary subsys- 
tem OA. This means that the probability V( m ) in each case 
(to) can be written as a product over the small sections of 
the boundary loops, therefore log 2 V( m ) is proportional to the 



length of the boundary. Since the combined boundary length 
of the cases (1) and (4) is equal to the combined boundary 
length of the cases (2) and (3), the topological Renyi entropy 



vanishes: 



C. Formula for the Renyi entropy 

Now we capitalize on the simplifications described above, 
and derive the Renyi entropy S^ A for an arbitrary state \ty) 
within the gauge theory. In the most general case, dA con- 
sists of n closed loops on the real lattice, and the loops have 
a combined length L. This means that they contain L spins 
and L stars acting on these spins (see Fig. [?}. Since there 
is a constraint on each loop due to the gauge structure, only 
L — n spins are independent. If we label these spins with 
1 < i < L — n, the 2 L ~" non-zero diagonal elements of pqa 
give the probabilities of | *]/) realizing the 2 L ~ n respective spin 
configurations {£f}. Since the projection operator onto the 
spin configuration {S| } is given by 2 n_i ]~J i (l + Sf of), the 
corresponding diagonal element reads 



(p 



1 



9A) SS 



2 L- 



L-r 



I + 



i=l 



(8) 



When expanding the product in Eq. (HJ and summing the 
squares of the resulting expressions for (pdA)w, the cross- 
terms cancel each other, and we obtain 



Tr [p 2 dA ] 



1 



2 L ~ 



{9i=0,l} 



n 



(9) 



where the sum is over all the 2 L ~™ configurations {^ = 0,1}, 
and hence over all possible products of the L — n independent 
spin operators of. If the edge occupied by the spin i con- 
nects the stars s^i and Sj^, the corresponding spin operator 
becomes af = A* x A x s . . In terms of the quasi-spin opera- 
tors A x s , the Renyi entropy then takes the form 



S dA = (L-n)-log 2 J2 ' 
{<?»=o,i} 



I*) 



(10) 



This expression has an entirely precise notation, but it is cum- 
bersome to use for calculating Sf A ■ To derive a more intuitive 
expression with a less precise notation, we expand the sum in 
Eq. ( [Tol l around the trivial configuration {qi =0}. Exploiting 
(A x ) 2 = 1, the Renyi entropy then becomes 



S. 



C 1 A = (L-n) -log 2 



£ <v&|i*if 2 i* 3 i* 4 |v]>) 



Sl,S2,S3,S4 



(ID 



where the sum inside the logarithm contains all 2 " possi- 
ble products with an even number of quasi-spin operators A x 
chosen from each closed loop of the subsystem dA. 
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To understand how Eq. (fTTT i works, we consider the two 
limiting ground states | f|~) and |0). In the first case, we have 

{it \A* 1 A%...A* a j 1\) = (t |<<x?...<x?| t) = i 

for all expectation values because | ff) has erf = +1 for 
all i. The sum inside the logarithm becomes 2 L ~", and the 
Renyi entropy 5f A vanishes, as expected for a product state. 
In the second case, |0) has A z s = +1 for all s, therefore 
(Q\A X S1 A X S2 . . .if 2r |0) = for all expectation values. The 
only exception is the trivial one: (0|0) = 1. The sum inside 
the logarithm is 1, and the Renyi entropy is S® A = L — n. 
When extracting the topological contribution, the terms oc L 
cancel because + L (4) = + (see Fig. 0}. On 
the other hand, the cases (1) and (4) have n = 2, while the 
cases (2) and (3) have n = 1, therefore the topological Renyi 
entropy is finite: S*J = 2. 

IV. PHASE TRANSITION IN THE QUASI-1D CASE 

In this section, we set k <C 1 in Eq. (HJ: this means that the 
spins on the vertical edges experience much smaller magnetic 
fields than those on the horizontal edges. The magnetic fields 
on the vertical edges only become important in the A ^> 1 
limit, therefore the phase transition between topological order 
and disorder at A ~ 1 occurs due to a competition between 
the star operators A z s and the magnetic fields on the horizontal 
edges. When investigating this phase transition, the magnetic 
fields on the vertical edges can be neglected (k — > 0), and the 
quasi-spin Hamiltonian in Eq. (O becomes 

H = -J2A*-\ £ A* S A%,. (12) 

s {s,s')eh 

Since there are only Ising couplings on the horizontal edges, 
this Hamiltonian is the direct sum of N independent ID TFIM 
copies along the horizontal chains of the lattice ll2lll . The sys- 
tem is therefore exactly solvable for all values of A. 

If we consider any of the independent ID horizontal chains 
and label the N stars on the chain with 1 < I < N, the Hamil- 
tonian of the corresponding ID TFIM reads 

N 

H = -J2{ A l + xA f A f-i)> (13) 
i=l 

where the periodic boundary conditions are taken into account 
by Ao = An. This Hamiltonian can be solved by following 
a standard procedure 1I22I . We first map the quasi-spins Ai to 
fermions via the Jordan- Wigner transformation 

if = 1 - 1cjc x , A- = f = eje™ c h , (14) 

where cj and q are standard fermionic creation and annihila- 
tion operators. The translational symmetry is then exploited 
by the Fourier transform 



where the sum is over the momenta k m = 7r(2m — i)/N with 
1 < m < N. We finally introduce new fermionic operators 
via the Bogoliubov transformation 

c fc = cos6»fc7fc +is'm8 k ~f j L k , (16) 
and the Hamiltonian in Eq. ( fT~3l > becomes 

F = ^A fe (2 7 | 7fe -l), (17) 

k 

where 7I and 7^ correspond to independent fermionic quasi- 
particles. The energies of these quasi-particles are propor- 
tional to Afc = J ef, + A 2 sin 2 k with e k = 1 — A cos fc, and 
the mixing angle appearing in the Bogoliubov transformation 
is 9 k = tan -1 [A sin k/ (e k + A fe )]. 

The ground state | O ( A)) of the Hamiltonian in Eq. ( TTZb is 
the direct product of N independent copies of the ID ground 
state |Oo). The ID ground state is defined by 7fc|S!o) = 
for all k, therefore its two-operator expectation values in the 
position representation are given by 

(c\c v ) = ^o|c i t C H^o> = ^E e ^ + ' ife '''^l c i^l^o> 

k,k' 

= ±Y e -*«- 1 ') sin 3 8 k , 

k 

( c '4> - ^E e " fe( ^' )cos2 ^' ( 18 > 

k 

(cicv) = ^5>-^- r >sin0 fc cos# fc , 

k 

(c\cl) = -±^ e -^-O sin ^ cos ^. 

To calculate the Renyi entropy, we need to evaluate the quasi- 
spin expectation values appearing in Eq. (fTTT i. These expecta- 
tion values are products of independent ID expectation values 
(flolA^A^ . . . Af 2r _ i Af 2 JQ ), where each pair of Af oper- 
ators can be expressed in terms of the fermionic operators as 

i'-i 

AfAf, = (cj + q) J] (l 2c] Cj ) (4 + c,) . (19) 

]=i 

Similarly, the quasi-spin expectation value in the Wilson loop 
for a D x D square region R becomes 

D D 

w R = (Ool Y[Af \n a ) D = (Ooin i 1 - 2c U) \n ) D - 

1=1 1=1 

(20) 

Using the identity 1 — 2cjq = (cj + q)(c| — q), the quasi- 
spin operator products appearing in both the Renyi entropy 
and the Wilson loop can then be written as simple products of 
cj ± q operators. On the other hand, the expectation values of 
these products can be reduced to the two-operator expectation 
values given in Eq. <fT~8T > by using Wick's theorem B23I1 . 
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FIG. 5: (Color online) Topological Renyi entropy S J as a function of 
magnetic field A below (a) and above (b) the critical point at Ac = 1. 
Five system sizes are plotted with d = D/3 for each: TV = 40 and 
D — 6 (i); TV = 80 and D — 12 (ii); TV = 120 and D — 18 (iii); 
TV = 160 and D = 24 (iv); TV = 200 and D = 30 (v). 
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FIG. 6: (Color online) (a) Reduced Wilson loop Wo as a function of 
region size D for A = 0.9 < Ac (red dashed line), A = 1.0 = Ac 
(black solid line), and A = 1.1 > Ac (green dotted line). The 
system size is TV = 500. (b) Reduced Wilson loop Wo as a function 
of magnetic field A for system sizes TV = 573 = 100 (dotted line) 
and TV = 5D = 500 (dashed line). The solid line is the result in Eq. 
J301 > for the thermodynamic limit. 



The exact dependence of the topological Renyi entropy on 
the magnetic field is plotted in Fig. [5] There are two phases: a 
topologically ordered phase at small A and a disordered phase 
at large A. These phases are separated by a clear phase tran- 
sition at A = Ac = 1, which coincides with the well-known 
critical point of the ID TFIM ll22tl . If we gradually increase A, 
the topological Renyi entropy drops to zero around Ac- This 
transition becomes sharper if we increase the system size TV 
as well as the dimensions D and d of the subsystems, there- 
fore we argue that 5j is discontinuous in the thermodynamic 
limit. The topological Renyi entropy is then constant in both 
limiting phases: the topologically ordered phase at A < Ac 
is characterized by Sj = 2, while the disordered phase at 
A > Ac is characterized by Sj = 0- 

The analogous exact behavior of the Wilson loop is illus- 
trated in Fig. [6] In the topologically ordered phase at A < Ac, 
the reduced Wilson loop Wo = approaches a finite 

constant in the D — > oo limit. This implies Wr cx exp(— j3D) 
and the presence of topological order. In the disordered phase 
at A > Ac, Wo decays exponentially with D. This implies 
Wr oc cxp(— f3D 2 ) and the absence of topological order. By 
looking at the dependence Wo (A) for a sufficiently large value 



of D, we can establish that the critical point separating the 
two different behaviors is indeed at Ac = 1. The results ob- 
tained for the topological Renyi entropy and the Wilson loop 
are therefore consistent with each other. 



V. PHASE TRANSITION IN THE ACTUAL 2D CASE 

In this section, we set k = 1 in Eq. <[TJ: this means that the 
spins on the horizontal and the vertical edges experience the 
same magnetic field. Up to an irrelevant additive constant, the 
quasi-spin Hamiltonian in Eq. (f3]i becomes 



H 



E 



i-A: 



A 



E 

(b,s>) 



A* A* 



(21) 



and the system is equivalent to the standard 2D TFIM. Since 
the Hamiltonian in Eq. ( f2Tb is not exactly solvable in general, 
we use perturbation theories around the exactly solvable lim- 
its at A = and A — > oo. The corresponding calculations are 
most efficiently performed by the method of perturbative con- 
tinuous unitary transformations (PCUT). The general method 
is discussed in the literature Ir20[ l24ll and we illustrate its use 
by the example of our particular problem. 



A. Perturbation theory at small magnetic field 

In the limit of A <C 1, it is useful to work in the quasi- 
spin representation because the unperturbed ground state |0) 
is then a product state. The perturbation theory is based on 
Eq. (1211 . where the second term is treated as a perturbation 
in the small parameter A « 1. Using the PCUT procedure 
described in Appendix [T] we obtain corrections to the Renyi 
entropy S^ A for each case of partitioning in Fig. [3] and the 
Wilson loop Wr for a square region R. The Renyi entropy 
after the first three corrections reads 



S. 



OA 



(L- 

K 
in2 



L 
In~2 



A^ 

4 



63A 4 503A 6 



27A 4 737A 6 



G4 



256 



04 



0(A 8 ) 



!)(> 



0(A 8 ) 



(22) 



where the boundary dA contains n closed loops with a com- 
bined length L and a total number of K corners that are suf- 
ficiently far away from each other. The analogous expression 
for the Wilson loop after the first three corrections is 



K 





[a 2 ; 


h 


J + ~ 


3A 4 


89A 6 


32 + 


128 + 



7697A 6 
3072 



0(A 8 



0(A 8 



(23) 



where the square region R has a boundary length L = AD and 
a corner number K = 4. 

Since K is merely an 0(1) constant, the corrections inside 
the exponential of Eq. ( 1231 are linearly proportional to the re- 
gion dimension D. The Wilson loop has therefore a functional 
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form Wr oc exp(— j3D) that shows the presence of topolog- 
ical order. Since the corrections to the Renyi entropy are all 
linearly proportional to either L or K, the corrections to the 
topological contribution oc n vanish. When calculating the 
topological Renyi entropy, the corrections cx L, K cancel be- 
cause the combined values of L and K in the cases (1) and 
(4) match those in the cases (2) and (3) (see Fig. [3). The 
topological Renyi entropy is therefore constant up to the third 
correction: — 2 + 0(A 8 ) in the A < 1 phase. 



B. Perturbation theory at large magnetic field 

In the limit of A ^> 1, it is useful to return to the physi- 
cal spin representation because the unperturbed ground state 
| ft ) is then a product state. Up to an irrelevant additive con- 
stant and an overall multiplicative factor A -1 , the 2D TFIM 
Hamiltonian in Eq. d2~TT i becomes 



(24) 



The perturbation theory is based on Eq. (1241 . where the sec- 
ond term is treated as a perturbation in the small parameter 
A -1 <C 1. Using the PCUT procedure described in Appendix 
[2l the Renyi entropy after the first three corrections is 



L 

ln~2 
K 
ln2 



A" 



A' 2 

32 + 1024 

35A- 6 



115A" 



4718592 



2359296 
0(A- 8 ) . 



0{\~ 



(25) 



Since the corrections to the Renyi entropy are all linearly 
proportional to either L or K, the corrections to the topo- 
logical contribution cx n vanish. The topological Renyi en- 
tropy is therefore constant zero up to the third correction: 
Sj = 0(A~ 8 ) in the A > 1 phase. 

To obtain a non-zero result for the Wilson loop expectation 
value (fi(A)|Wfl|fi(A)), we need to consider higher orders of 
perturbation theory. Since Wr is a product of D 2 star op- 
erators Ag , the first non-zero contribution to Wr appears at 
order D 2 /2 in perturbation theory. At this order, Wr links 
order D 2 /2 states to each other, therefore Wr oc A _d . This 
result can be rearranged into the form Wr oc exp[— ln(A) D 2 ] 
that shows the absence of topological order. 



C. Discussion of the phase transition 

The results of the perturbation theories indicate two distinct 
phases around the limits A = and A — > oo. The phase at 
A <C 1 is topologically ordered because the topological Renyi 
entropy is non-zero and the Wilson loop follows a perimeter 
law: Wr oc exp(— (3D). Conversely, the phase at A S> 1 is 
disordered because the topological Renyi entropy is zero and 
the Wilson loop follows an area law: Wr oc exp(— (3D 2 ). 

The topological distinctness implies at least one phase tran- 
sition between the two limiting phases, and we argue that there 



can only be one phase transition. Recall that the TCM with ex- 
ternal field is equivalent to the standard 2D TFIM when k — 1 . 
In particular, the quantities Sj and Wr that describe topolog- 
ical order can be expressed in terms of the 2D TFIM correla- 
tion functions. A phase transition is therefore only possible at 
the critical point of the 2D TFIM, which has been determined 
by various numerical methods ll25ll to be at Ac ~ 0.33. This 
critical field is also consistent with previous numerical studies 
on the TCM with external field [fll [ll . 

It is clear that the perturbation theories around the two lim- 
its need to break down at A = Ac- On the other hand, the re- 
sults of the perturbation theories hold because the expansions 
in Eqs. 1221 . d23"T l. and d25l l have particular structures: they 
each contain two power series in A that are proportional to the 
boundary length L and the corner number K. It is plausible 
that higher order corrections preserve this form and only add 
further terms to the respective power series. Terms that are 
not linearly proportional to either L or K only appear when 
the order of the perturbation theory exceeds the dimensions 
D and d of the subsystems (regions). Since these dimensions 
are macroscopic in the thermodynamic limit, the perturbation 
theories can only break down at infinitely large orders. These 
in turn become important at the radii of convergence where 
the series actually diverge. If we write the power series in 
Eqs. d22l and 1231 as Y^kLi a kX 2k and those in Eq. l25l > as 
SfeLi bkX~ 2k , the critical field Ac is given by 



lim 



a k 



ak+i 



lim 



3fc+l 



bk 



(26) 



Although it is not possible to determine these limits from a 
finite-order perturbation theory, we can give estimates for the 
critical field by looking at the first couple of terms and calcu- 
lating analogous quantities. The resulting estimates are sum- 
marized in TablelU they suggest 0.2 < Ac < 0.5. This range 
is fully consistent with Ac ~ 0.33. 



Estimates for Ac 




\J 02/03 


Eq. 422j 


Series oc L 
Series oc K 


0.504 


0.433 
0.383 


Eq. ([23) 


Series oc L 
Series oc K 


0.500 


0.447 
0.367 




Estimates for Ac 


^62/61 


\Jb-ijb2 


Eq. (25) 


Series oc L 
Series oc K 


0.177 


0.223 



TABLE I: Estimates for the critical field Ac obtained from the power 
series of Eqs. J22K l !23t . and d25t . 

The most remarkable result of this section is that the topo- 
logical Renyi entropy is constant in both limiting phases: 
S 2 = 2 in the topologically ordered phase and S J = in 
the disordered phase. This happens because the perturbative 
corrections to do not contain any topological contribu- 
tions oc n in Eqs. d22l and (f25j. The topological Renyi en- 
tropy is therefore an exclusive function of the phase: it can 
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only change if a phase transition takes place. We argue that 
S2 is a good probe of topological order with the potential to 
characterize topologically ordered phases. 



VI. COMPARISON AND DISCUSSION 

A. Perturbation theories in the quasi-lD case 

Although the quasi- ID case is exactly solvable, it is instruc- 
tive to treat it with perturbation theories as well: the results 
obtained this way are directly comparable with those in the 
actual 2D case. Using a modified version of the PCUT pro- 
cedures described in the Appendix, we find analogous expres- 
sions to those in Eqs. ( l22b . d23l ), and d25l l. Without includ- 
ing the detailed calculations, the Renyi entropies after the first 
three corrections in the two limiting regimes are 



A^ TA^ 5A^ ^/»s\ 

1 h h O A 8 

4 64 96 1 ' 



H' 
ln~2 



5A^ 3A^ _ . . . 
^4- + ^ + 0(A) 



Estimates for Ac 


\/ai/a 2 




E q .m Series « L ' 

Series oc H' 


1.512 


1.449 
0.913 


Eq. (|29l Series oc 1 


1.414 


1.225 




Estimates for Ac 


y/b 2 /bi 




Series oc H' 


0.500 
0.661 


0.700 
0.798 



(A « 1), (27) 



TABLE II: Estimates for the critical field Ac obtained from the 
power series of Eqs. ( 127) . ( 128) . and ( 129) . 



Ac marks the breakdown of the perturbation theories, and it is 
again related to the appropriate radii of convergence. The es- 
timates obtained with the method of Sec. IV Cl are summarized 
in Table HI1 they suggest 0.5 < Ac < 1-5. This range is fully 
consistent with Ac = 1. 



B. Comparison with the actual 2D case 



V 
hi~2 
H' 
In~2 



A^ 

8 
A' 2 
8 



A" 4 

7A- 4 
128 



47A~ 6 
3072 
107A- 6 



3072 



+ 0(A" 8 ) 
0(A" 8 ) 



(A » 1), 



(28) 



where L' is the combined horizontal length of the boundary 
dA, and H 1 is the number of horizontal sections with a non- 
zero length contributing to V . Since the corrections to the 
Renyi entropies are all linearly proportional to either L' or 
H', there are no topological corrections oc n. The topological 
Renyi entropy is therefore constant Sj = 2 + 0(A 8 ) at A <C 1 
and constant Sj = 0(A~ 8 ) at A > 1. 

When taking into account the first three corrections, the re- 
duced Wilson loop in the A <C 1 regime becomes 



Wo = exp 



(29) 



which indicates Wr oc exp(— j3D) and the presence of topo- 
logical order. In the A ^> 1 regime, the first non-zero con- 
tribution to Wq appears at order D/2 in perturbation theory. 
This contribution is Wq oc A~ d = exp[— ln(A) D], which 
indicates Wr oc exp(— (3D 2 ) and the absence of topological 
order. Note that the power series inside the exponential of Eq. 
(|29l suggests that Wo takes the exact form 



Wn 



exp 



E 



2fc 



A 

Ik 



A 2 ) 1/4 



(30) 



in the thermodynamic limit. This result is consistent with the 
critical field Ac = 1 obtained from the exact treatment. 

The perturbative expansions in Eqs. (|27j , (|28l , and ((29) 
each contain at least one power series in A. The critical field 



When discussing the phase transition in the actual 2D case, 
we argued that it occurs at the critical point Ac ~ 0.33 of 
the equivalent 2D TFIM and that the two limiting phases are 
characterized by different constant values of the topological 
Renyi entropy. The argument only referred to the perturbation 
theories and the equivalence with the 2D TFIM. On the other 
hand, the quasi- ID case is more versatile because an exact 
solution is available. The exact treatment of the quasi-lD case 
suggests a behavior that is entirely analogous to our claims for 
the actual 2D case: the phase transition occurs at the critical 
point Ac = 1 of the equivalent ID TFIM, and the topological 
Renyi entropy is constant in the two limiting phases. 

A direct comparison between the respective perturbation 
theories also provides evidence that the ID and the 2D sys- 
tems are similar in terms of their phase transitions. The be- 
haviors of the Ac estimates and their relations to the actual Ac 
are entirely analogous in the two cases. First, the estimates 
are all reasonably close to the actual Ac- Second, the esti- 
mates converge towards Ac as the order is increased. Third, 
the estimates from the A « 1 series generally overestimate, 
while those from the A ^> 1 series underestimate Ac- These 
similarities suggest that the phase transitions in the ID and the 
2D cases are analogous, therefore the conclusions drawn from 
the exact treatment in the quasi-lD case are applicable to the 
physically more interesting actual 2D case as well. 



VII. CONCLUSIONS 

In this paper, we investigated the quantum phase transition 
between the topologically ordered and the disordered phases 
of the TCM with external magnetic field. The variation in 
topological order was probed via Sj: the topological Renyi 
entropy of order 2. We determined the exact field dependence 
of S2 in the computationally simpler case (quasi- ID case) 
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and established perturbation theories in the physically more 
interesting case (actual 2D case). It was demonstrated that Sj 
takes distinct values in the two phases and has a discontinu- 
ity at the quantum phase transition. We therefore argue that 
S2 is a good probe of topological order that can effectively 
characterize topologically ordered phases. 

The equivalence between the quasi- ID case of our prob- 
lem and the exactly solvable ID TFIM is a quite remarkable 
tool for obtaining exact results. So far it has provided us with 
an exact treatment of the quasi- ID case and a corresponding 
exact <5j(A) dependence. In perspective, such an exact treat- 
ment also makes it possible to search for critical exponents 
that can reveal the topological character of the quantum phase 
transition. Moreover, the exact time dependence of the system 
far away from equilibrium can be studied, as for example, in 
the case of a quantum quench lf26l l27ll . 

It is important to point out that the Hamiltonian in Eq. (Q3 
preserves the Z2 gauge structure of the bare TCM for all val- 
ues of the magnetic field A. This gauge structure justifies the 
simplifying step of substituting the subsystem A by its bound- 
ary dA when calculating Sj (thin subsystem). Indeed, as 
long as the gauge structure is preserved by the perturbation, 
the ground state can be expressed as a superposition of loop 
configurations. For such a system, all the relevant topological 
constraints are necessarily connected to the subsystem bound- 
ary dA. For example, in our Z2 gauge theory, the topological 
constraint manifests itself in the fact that there are an even 
number of spins flipped on each boundary loop of dA. On 
the other hand, considering only the boundary is the essential 
simplification we need for deriving the Renyi entropy formula 
in Eq. (fTTT l. which in turn makes the exact treatment in the 
quasi- ID case and the perturbation theories in the actual 2D 
case possible. The gauge structure also explains why the topo- 
logical Renyi entropy is conserved du ring a quantum quench 
with a gauge-preserving Hamiltonian II27II . 

For a more generic Hamiltonian, the Z2 gauge structure 
is broken. This means that the spin configurations with an 
odd number of spins flipped on a boundary loop of dA are 
allowed, therefore the topological constraint is no longer 
connected to the subsystem boundary dA. Note that we 
can also achieve an effective gauge structure breaking by 
drawing the boundary loops of dA on the dual lattice rather 
than on the real lattice (see Fig. |4]i because they can then 
intersect with the spin loops on the dual lattice at an arbitrary 
number of points. To recover the robustness of S 2 in such a 
non-gauge-preserving case, one needs to calculate it by using 
the original subsystem A (thick subsystem). This complicates 
the situation because the reduced density matrix pA is not 
diagonal and so Eq. (fTTT l becomes invalid. However, we 
believe that if a generalization of the Renyi entropy formula 
is found, the results in this paper can be extended to the more 
generic non-gauge-preserving case as well. This further step 
is crucial for verifying the robustness of the topological Renyi 
entropy against generic perturbations and hence proving its 
applicability as a non-local order parameter for topologically 
ordered phases. 
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Appendix: Detailed descriptions of the PCUT calculations in the 
actual 2D case 

1. PCUT calculation at small magnetic field 

When considering Eq. (fZTT i in the A <C 1 limit, we can use 
the PCUT procedure to relate the eigenstates of the perturbed 
Hamiltonian H with A > to those of the unperturbed Hamil- 
tonian Ho with A = 0. This method relies on the concept of 
elementary excitations. In the case of Hq, these excitations 
are flips of stars (quasi-spins) A z s with an energy cost of 2 for 
each. They appear pairwise when switching on the perturba- 
tion, and the perturbed Hamiltonian can be written as 

H = 2Q + f +2 +f + T_ 2 , (A.l) 

where Q counts the number of excitations, and T n is the com- 
ponent of the perturbation that changes the number of exci- 
tations by n. It can be verified that [Q,T n ] = nT n and that 
T,! = T_ n . The explicit forms of the terms in Eq. (IA.1) are 

Q = \ E - A s) > f +2 = -A E A-A-„ (A.2) 

8 (S,S>) 

To = -A J2 (AtAj, +AjA+) , f_ a = -Xj2 A tAf„ 

(s,s>) (s,s>) 

where Af = (A*±iA*)/2 are the standard spin raising and 
lowering operators. In the basis of the Hq excitations, the term 
Q is diagonal, while the terms T n are non-diagonal. The ap- 
plication of the PCUT involves an iterative sequence of steps 
to construct a unitary basis transformation U(l) such that the 
transformed Hamiltonian H(l) = W (l)HU(l) changes con- 
tinuously from H at I = to a block-diagonal form at I 00. 
The blocks in the asymptotic form H' = H(oo) correspond to 
subspaces of constant excitation number, and the excitations 
can be found by solving the blocks. Note that these excitations 
belong to the perturbed Hamiltonian H, therefore they are not 
the same as the original Hq excitations. To avoid confusion, 
we refer to them as quasi-excitations. 

According to the standard procedure of the PCUT 1I20II24I1 . 
we write H{1) = 2Q + f +2 (l) + f (l) + T_ 2 (0 as in Eq. 
(lA.ll i, and define fj(l) = T +2 (l) - T_ 2 (Z). If we then require 
U(l) to satisfy the equation d t U(l) = -U(l)fj(l), it follows 
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that<9/H(Z) = [f)(l),H(l)]. In terms of the components T n (l), 
this equation for H(l) becomes 



diT (l) = 2 [T +2 (0,T_ 2 (0 
dtf +2 (!) = -4T +2 (l)+[f +2 (l),T (l) 
dif_ 2 (l) = -4T_ 2 (/)+ [f o (0,r_ 2 (0 



(A.3) 



The last two equations show that T +2 (oo) = T_ 2 (oo) = 0, 
which is consistent with the block-diagonal form of H(oo). 
To solve the equations for U(l) and T n (l) iteratively, we write 
these quantities in a series as 

oo oo 

U(1) = J2UW(1), f n (l) = J2Ti k \l). (A.4) 



k=0 



k=l 



The equations for U (I) and T n (l) then take the forms 



k-l 



8^(1) 



j2^ j Hi){fi k 2 - j \i)-f^- j \i)}, 



3=0 
k-l 



2£ [f« (0,^(0 
-4fW(0 + 2[fS(0,r o ^(0 
-4fW(0 + ^[fW(0,f^(0 



(A.5) 



and the corresponding starting conditions at I = become 

1 (k = 0) 



w I (fc > 1) 



(A.6) 



T«(0) = { 



T, 
o 



(fc = i) 

(jfc > 2) 



If we apply the PCUT up to second order in A, the relevant 
terms in the series of Eq. (I A. 4) are 



^(0 



i, ^(0 = ^(T +2 -f_ 2 )(i- e - 4 '), 

:i- e - 4 *) 2 (A?) 



32 
1 

16 



T +2 - T_; 



T+ 2 + T_ 2 , To 



1 - (1 + 40 e 



T (1 ^) = To, T,T(0 



(2) 



r +2 ,T_, 



l-e 



-41 



-8/\ 



f^(0 = T ±2 e- 4 ', 2^(0= ± 



(2), 



T± 2 ,T 



the basis transformation J7 = C/(oo) becomes 



1 

Tg 



To,T_ 2 



T+2, To 



32 



+2-L+2 + 1-21-2 — 1+2I-2 — l-2l- 



21-2 — 1+21-2 



-21+2 , 



(A.8) 



H' = 2Q + f 



T +2 ,T_ 5 



(A.9) 



and the asymptotic Hamiltonian takes the form 

1 r. 
4 

which indeed conserves the number of quasi-excitations. The 
same procedure can be continued to arbitrary order in A, but 
the calculations quickly become cumbersome. 

Since the ground state |f2(A)) of H is the only state with no 
quasi-excitations, it has its own block in H' . To express this 
state in terms of the physically transparent Hq excitations, we 
use the basis transformation: |0(A)) = U\Q). When calculat- 
ing the first two perturbative corrections to the ground state at 
A = 0, the perturbed state |0(A)) needs to be properly nor- 
malized up to A 4 . Applying the PCUT up to fourth order with 
the aid of a computer, the perturbed ground state becomes 



MX)) = u\o) 

A 

4 ~ 



2 \ 2 



1 



(A 2 



N 2 X 
16 
15)A 3 



(A 4 - 95A 2 )A 4 
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G4 
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X X 
X X 



Eh 

2N 2 

^E 



(A. 10) 



A 2 
16 



2N 2 

E 

2 JV 4 -9JV 2 



where the equivalent states related to each other by transla- 
tional and rotational symmetries are labeled by the relative 
positions of the star excitations (x) in them, and the number 
of states in each equivalence class is given by the number be- 
low the corresponding sum. The notation [. ..][.. .] means that 
there are two clusters of excitations that are independent of 
each other: they are not in a relative position characterizing 
any other equivalence class. 

The ground state in Eq. dA.lOl i is indeed properly normal- 
ized up to fourth order in A because 



(0(A)|0(A)> 



1 



2A 2 



2N 2 



TV 2 A 2 (A 4 - 95A 2 )A 4 



16 512 

A _ (A 2 - 15)A 3 n2 
4 

A- 



G4 



(A.ll) 



+ (2A 2 + A 2 ) (- 



+ (2A 4 -9A 2 ) f^J =1 + 0(A 6 ). 

To calculate the Renyi entropy, we need to evaluate the ex- 
pectation values of the products appearing in Eq. (fTTl i for the 
ground state. The expectation values having a contribution up 
to A 4 to the Renyi entropy are 
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(A. 12) 

f 0(A 4 ), 

f 0(A 4 ), 
0(A 4 ), 
0(A 4 ), 



where the notation is analogous to that in Eq. ( lA.lOl i. For 
example, ( x X } = (f2(A)|Ag |0(A)), where s and s' are 
any two nearest-neighbor stars. If the boundary dA consists of 
?i closed loops with a combined length L and a total number 
of K corners that are sufficiently far away from each other, 
the Renyi entropy is given by 
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(A.13) 



The perturbative corrections are linearly proportional to either 
L or K, and they are independent of n. 

Now we consider a Wilson loop for a square region R with 
boundary length L = AD and corner number K = 4. Accord- 
ing to Eq. the Wilson loop is Wr = 1 for the unperturbed 
ground state |0) because A z s = +1 for all s. The structure of 
the perturbed ground state |f2(A)) in Eq. (IA. 10b shows that 
Wr can only be —1 instead of +1 if an odd number of excita- 
tions are inside R. Taking into account all possibilities up to 
fourth order in A, the Wilson loop becomes 
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FIG. 7: Equivalence classes for the relative positions of the star ex- 
citations (x) at the level of the third corrections (at order A 6 ). 



The expressions in Eqs. JA. 13t and dA.14b give the first two 
corrections to A and Wr. With the aid of a computer, the 
third corrections proportional to A 6 can be found in a similar 
manner. In this case, the state |0(A)) must be properly nor- 
malized up to sixth order in A, and one needs to consider all 
the relative excitation positions shown in Fig. [7] Without in- 
cluding the detailed calculations, the final results for the Renyi 
entropy and the Wilson loop are 
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(A. 14) 



The perturbative corrections inside the exponential are lin- 
early proportional to either L or K. 



The features noticed after the first two corrections remain in- 
tact after the third corrections as well. 



2. PCUT calculation at large magnetic field 

When considering Eq. d24l > in the /i e A -1 « 1 limit, the 
PCUT procedure is entirely analogous to the one described in 
Appendix [T| The elementary excitations of the unperturbed 
Hamiltonian with fi = are flips of physical spins of with 
an energy cost of 2 for each. The perturbed Hamiltonian with 
fx > can be written as 

H = 2Q + f +4 +f +2 + f + f- 2 +T- 4 , (A.17) 

where the respective terms take the explicit forms 

Q = ^Et 1 -^)' ^ A - 18 > 

i 

s ± igs 
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The sum in ± contains all inequivalent products of the four af 
operators in which the number of the af factors is 2 — n/2 
and that of the af factors is 2+n/2. Applying the PCUT with 

fj(l) = f +4 (l) + f +2 (l) - f_ 2 (0 - T- 4 (0 up to fourth or- 
der, we find that the perturbed ground state which is properly 
normalized up to p 4 is given by 



|Q(A)) 



2 ,,2 



1 



N Z p 
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62N< 



4 1 



32768 



r, 



a _ { N i _ 1 \ 

3 



1024 



E 

N 2 




48 



2N 2 



(A.19) 



64 



E 

i(JV 4 -5JV 2 ) 



where the equivalent states are labeled by the relative posi- 
tions of the spin excitations (o) in them [cf. Eq. (IA.10H . 

We now consider a subsystem with a boundary dA con- 
taining L spins and L stars acting on these spins. The total 
number of corners is K as above. The diagonal elements of 
the density matrix pqa can be obtained directly in the basis of 
the physical spins af. The element corresponding to af = +1 
for all L spins (no star excitations on the boundary) is 



Note that there are L ways of choosing two neighboring spins 
from dA. Since the contribution of the remaining diagonal 
elements is 0(p 6 ) to the Renyi entropy, we find that 



Tr [p 2 9A ] 



{P8a)1q + L (pdA)li 



L , / L 

l -T2» 



V1024 2048 
and the Renyi entropy takes the form 



S d 2 A 



1 

In~2 



— p 2 + —p 4 + 0(p 6 ) 
32 P 1024 



(A.22) 

0(/i 6 ), 

(A.23) 



The perturbative corrections are again linearly proportional to 
the boundary length L. Furthermore, there are no terms oc K 
in the first two corrections determined here. 

With the aid of a computer, the third correction to S® A can 
be calculated in a similar manner. In this case, the state |f2(A)) 
must be properly normalized up to sixth order in p, and one 
needs to consider all the relative excitation positions shown in 
Fig. [8] Without including the detailed calculations, the final 
result for the Renyi entropy is 



5. 



OA 



L 

h^2 
K 
In~2 



32 + 1024 
35A" 6 



115A- 
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2359296 
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(A.24) 
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(A.20) 
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N 2 (2L + 
2 ( 5i 



(L 



V 8192 
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p 4 + 0{p 6 



while the one corresponding to af = — 1 for any two neigh- 
boring spins and af = +1 for the remaining L — 2 spins (one 
star excitation on the boundary) is 



(a 



dA, 
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EL 

64 



0(p 4 



(A.21) 



After the third correction, there is a corner contribution cx K, 
but still no topological contribution oc n. 
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FIG. 8: Equivalence classes for the relative positions of the spin ex- 
citations (o) at the level of the third corrections (at order p 6 ). 
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